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Abstract 

In this paper, we get the ID approximate analytical solution of the plasma electrostatic 
wake driven by the laser, and get the oscillating frequency of the wake modified by the 
nonlinear high order terms, and find that the frequency depends on the oscillating amplitudes 
p— I of the wake and the laser which is the general nonlinear phenomenon. Finally we analyze the 

longitudinal beam dynamics in this electrostatic wake, and find that the high order terms 
II don't change the topology of the longitudinal phase space. 
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1 Introduction 

O 

• ^ Laser plasma wakefield accelerators (LWFA) , in which the plasma wakefield is excited by an 

intense laser to accelerate the particles, have demonstrated accelerating gradient of hundreds of 
GV/m, and have become the highlights of advanced accelerator concept [1]. In this paper, we 

L^j* get the ID analytical solution of the plasma electrostatic wake driven by the laser, and get the 

modified oscillating frequency of this wake. Finally we analyze the longitudinal beam dynamics 
^ in this electrostatic wake, and find that the high order terms don't change the topology of the 

CN longitudinal phase space. 

On 
^ — I 

^ 2 Basic equations for LWFA 

cn 

We consider the homogeneous unmagnetized plasma. Laser travels in the z direction with 
^ the initial frequency wq f^p i is the electron plasma frequency, and the group and the phase 

'y^ speed Vg ^ Vp ^ c {c is the light velocity ) . We adopt one dimensional theory, the LWFA system 

is governed by the following set of equations [2] 
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where, E, B are the electromagnetic fields, and v is the electron velocity. We introduce the 
coordinates ^ = kp{z — ct){kp = tOp/c), and normalized scalar and vector potentials (/)(^) and a{^) 
. Such that E'li = — = ^^§| , where _L and || refer to the components perpendicular 

and parallel to z . From the equations above, we can get the equations about (p{S,) and a(^) , 
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where kp = is the electron plasma wave number, 7 is the relativistic factor of electron. From 
the Hamiltonian of the electron and the continuity equation, one can get the following relation 

ne ^ _J_ ^ l + a^ + (l + 0)^ , . 
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Finally, we get the ID self-consistent equations for LWFA system 
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the difference between the equations above and that in [2j is that we use the direct differential 
and the equations are self-consistent, we don't need other variables to solve the equations. 



3 The oscillating frequency of the wake 

In this paper, we assume that a is unchanged and (j) is in the vicinity of zero (small variable) 
</> oc e to get the analytical solution of the electrostatic plasma wake. We expand the equation 
(9) into the form 

(/>" + (/>- + 203 _ ^ ^ ^2^2^-^ _ 2^ + 302 _ ^ 5^4 ^ )^ (^g) 

Using KBM method [3j, one expand as following 

N 

= Mcos0 + ^e>„(n,0) + O(e"). (11) 
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The equations of u and 6 are 

^ = j;6M„(n), (12) 
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_ = ^e"6"e„(n). (13) 

From the equation (10) , we remove the secular terms to obtain the first order solution 

Ai = 0,6*1 = (14) 

37/2 2 

= - \cos2e, (15) 

and also the second order solution 

^2 = O,02 = a2-^, (16) 

<p2 = a^ + ^cos30. (17) 
So the solution of the electrostatic plasma wake potential to second order is 

(j) = ucosO + e(— -cos2e) + e'^{a'^ + ——cosSe), (18) 

4 4 128 

where 6 = kp{l + - ^){z - ct). In (18), if we reserve the zero order term only, then , the 
wake oscillating frequency is the usual electron plasma frequency ujp. But the wake oscillating 
frequency is modified as (19), if we include the high order terms, the modified frequency uj'p is 
related to the oscillating amplitude 

Q7/2 

^; = Ml + «'-^)- (19) 

The modified frequency depends on the square of the amplitude, this is the general nonlinear 
phenomenon |14j and can be extended to the more common case in plasma physics. 

4 The longitudinal beam dynamics 

Now we begin to analyze the longitudinal beam dynamics in the electrostatic potential . In 
doing so, the expression of should be expressed as = kp{l + — ^)(-z — Vpt),Vp < c. We 
incorporate the new variable s = kpZ, the longitudinal equations of motion are 

where (3p = Vp/c, /? = (1 — 7^^)2 . The function cj) is shown in the figure 1, 
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Figure 1. The function (f) ( a = 0.5 ). 
The fixed points of the longitudinal equations of motion are 7 = 7p, = nvr ( n is integer 

1 



From the Jacobian matrix of (20) and (21), the characteristic value A satisfies 



It 



is obvious that = is the stable point and 6 = — vr, vr are the saddle points. If we remove the 
second and third terms in (18), we will get the same conclusion. So the high order terms in the 
electrostatic potential don't change the topology of the longitudinal phase space. Figure 2 gives 
the longitudinal phase space near the zero phase, when u = 0.5, f3p = 0.98. 
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Figure 2. The longitudinal phase space near zero. 
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Figure 3. The whole longitudinal phase space. 
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